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We apply the effective field theory approach to the three-nucleon system. In partic- 
ular, we consider neutron-deuteron scattering and the triton. Precise predictions 
for S = 3/2 scattering are obtained in a straightforward way. In the S = 1/2 
channel, however, a unique nonperturbative renormalization takes place which re- 
^ ' quires the introduction of a three-body force at leading order. We also show that 

VP invariance under the renormalization group explains some universal features of the 

, three-nucleon system. 

o 
in , 

, 1 Introduction 

0^ ' Effective field theories (EFT) are a powerful concept designed to explore a 

separation of scales in physical systemsu For example if the momenta k of 
two particles are much smaller than the inverse range of their interaction 1 / R, 
observables can be expanded in powers of kR. Following the eaxJy work of 
Weinbergp EFT's have become quite popular in nuclear physics Era Their ap- 
plication, however, is complicated by the presence of shallow (quasi) bound 
states, which create a large scattering length a ^> R. In this finely tuned 
case, the perturbative expansion in kR already breaks down at rather small 
momenta k ~ 1/a. In order to describe bound states with typical momenta 
k 1/a, the range of the EFT has to be extended. A certain class of diagrams 
has to be resummed which generates a new expansion in kR where powers of 
ka are kept to all orders. For the two-nucleon system^Liiower counting that 
incorporates this resummation has been found recently&ErQ Pions are included 
and treated perturbatively in this scheme. It has successfully been applied to 
NN scattering and deuteron physics □ 

The three-nucleon system is a natural testground for the understanding of 
the nuclear forces that has been reached in the two-nucleon system. However, 
the extension of these ideas is not straightforward, as the three-nucleon system 
shows some remarkable universal features. It has been found that different 
models of the two-nucleon interaction that are fitted to the same low-energy 
two-nucleon (NN) data predict different but correlated values of the triton 
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binding energy B3 and the iS-wave nucleon-deuteron (Nd) scattering length 
ag 1 ^ in the spin 5=1/2 channel; all models fall on a line in the B3 x a^ 2 ^ 
plane, the Phillips lineEl Other universal features of three-body systems are the 
existence of a logarithmic spectrum of bound states that accumulates at zero 
energy as the two-particle scattering length 02 increases (the Efimov effect El), 
and the collapse of the deepest bound state when the range of the two-body 
interaction T2 goes to zero (the Thomas effectllil). 

We will show that the renormalization of the three-nucleon system is non- 
perturbative and requires a three-body force at leading order. In addition 
to the description of experimental data from a few low-energy parameters, 
to understand the above mentioned universal features in a unified 



2 Neutron-Deuteron System 

In order to avoid the difficulties due to the long range Coulomb force, we con- 
centrate here on the neutron-deuteron (nd) system. We also restrict ourselves 
to scattering below the deuteron breakup threshold where S-waves are domi- 
nant. There is only one scale at low energy, \JMB2 ~ 40 MeV, where B2 is 
the binding energy of the deuteron and M the nucleon mass. Since yJMBi is 
small compared to the pion mass, the pions can be integrated out and only the 
nucleon degrees of freedom remain. The leading order of this pipnless theory 
is equivalent to the leading order in the KSW counting scheme B where pions 
are included perturbatively. As a consequence, the extension of our results to 
higher energies is well defined. Recently, a number of two nucleon observables 
have been studied in the pionless theory as wellE-3 

There are two S'-wave channels for neutron-deuteron scattering, corre- 
sponding to total spin S = 3/2 and S — 1/2. For scattering in the S — 3/2 
channel all spins are aligned and the two-nucleon interactions are only in the 
3 Si partial wave. The interaction is repulsive and the Pauli principle forbids 
the three nucleons to be at the same point in space. As a consequence, this 
channel is insensitive to short distance physics and very precise predictions are 
obtained in a straightforward way. There is also no three-body bound state 
in this channel. The S — 1/2 channel is more complicated. The two-nucleon 
interaction can take place either in the 3 S\ or in the 1 So partial waves. This 
leads to an attractive interaction which sustains a three-body bound state, the 
triton. The S — 1/2 channel also shows a strong sensitivity to short distance 
physics as the. Pauli principle does not apply. As a consequence, it displays 
the Thomas O and Efimov li!J effects. The generic features of this channel are 
very similar to the system of three spinless bosons. There is a strong cutoff 
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dependence even though all Feynman diagrams are finite. As we will shflaz- 
the renormalization requires a leading order three-body force countertermlialiJ 
Since the details oLthe renormalization in the three-body system are discussed 
in Bedaque's talkf^l we will mainly focus on the application of these ideas to 
the nd system. 

Let us start from the assumption that the three-body force is of natural 
size. The lowest order effective Lagrangian is then given by 



C = N ] (id + ^— \ N -Cq (N T T 2 aa 2 N) ] ■ (N T T 2 aa 2 N) (1) 




-C S (N T o 2 tt 2 N^ ■ (N T a 2 TT 2 N) + ... , 

where the dots represent higher order terms supressed by derivatives and more 
nucleon fields, a (r) are Pauli matrices operating in spin (isospin) space, 
respectively. The contact terms proportional to Cq (Cq) correspond to two- 
nucleon interactions in the 3 Si ( 1 5'o) NN channels. Their renormalized values 
are related to the corresponding two-body scattering lengths a\ and a 2 by 
Cq' = Tra,2 /(2M). Since no derivative interactions are included, this La- 
grangian generates only two-nucleon interactions of zero range. For practical 
purposes, it is convenient to rewrite this theory dry introducing "dibaryon" 
fields with the quantum numbers of two nucleonsE£l In our case, we need two 
dibaryon fields: (i) a field T with spin (isospin) 1 (0) representing two nucle- 
ons interacting in the 3 <Si channel (the deuteron) and (ii) a field S with spin 
(isospin) (1) representing two nucleons interacting in the channel. Using 
a Gaussian path integration, it is straightforward to show that the Lagrangian 
from Eq. ([!]) is equivalent to 

C = fjd + J^JJV + Arft-f + A s S f -S (2) 

-— (ft ■ N T T 2 aa 2 N + h.c^j - — (V ■ N T g 2 tt 2 N + h.c. 

At first it may look like the Lagrangian, Eq. (^), contains more parameters 
than the original one, Eq. (|l|). However, the scales At and As are arbitrary 
and included in Eq. (^) only to give the dibaryon fields the usual mass dimen- 
sion of a heavy field. They can easily be removed by rescaling the dibaryon 
fields. All observables depend only on the ratios g\ s /At,s- Since the theory 
is nonrelativistic, all particles propagate forward in time, the nucleon tadpoles 
vanish, and the propagator for the nucleon fields is 

i 

Po - p /2M + le 



3 



The dibaryon propagators are more complicated because of the coupling to two- 
nucleon states. The bare dibaryon propagator is simply a constant, i/ As,t, 
but the full propagator gets dressed by nucleon loops to all orders as illustrated 
in Fig. [I]. The nucleon loop integral has a linear UV divergence which can 




Figure 1: Dressing of the bare dibaryon propagator. 

be absorbed in g 2 -, s j A^s, a finite piece determined by the unitarity cut, and 
subleading terms that have already been omitted in Eq. (^). Summing the 
resulting geometric series leads to 

where <7t,s & n d Ay 5 now denote the renormalizcd parameters. The NN scat- 
tering amplitude in the respective channel is obtained by attaching external 
nucleon lines to the dressed propagator. In the center of mass frame, the 
S-wave amplitude in the 3 Si, 1 So channels for energy E = k 2 /M is 

3 ' 1 T(k) = ^[-=^-ik\ , (5) 
and the renormalized parameters gx,s and At,s can be determined from 

s.t M 9l,T (f ,s 

3 S = 3/2 rid-Scattering 

For S = 3/2 scattering only the dibaryon field T contributes. The first few 
diagrams that contribute are shown in the first line of Fig. ^. Since the 
leading piece of all diagrams is of the order ~ Mg^/Q 2 , they have to be 
resummed for typical momenta Q ~ X/a\- This is conveniently achieved by 
solving the integral equation given by the second line in Fig. || Performing the 
integration over the time component of the four-momentum flowing in the loop 
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Figure 2: The first few diagrams contributing to S = 3/2 nd-scattering (first line) and the 
integral equation which sums the diagrams to all orders (second line). 



and projecting onto the S-waves, we obtain the scattering amplitude a(fc,p)0 



- A (l/a* + V3p 2 /4 - ME) 1 a{p,k) 
q 2 dq 



(7) 



-K{p, k) 



q 2 — k 2 — ie 



K(p, q) a{q, k) , 



where 



2pq 



q +pq- p 



ME 



pq — p* 



ME 



(8) 



and ME = 3fc 2 /4 — l/( a 2) 2 ^ s the total energy. k(p) denote the incoming 
(outgoing) momenta in the center of mass frame. The amplitude a(p, k) is 
normalized such that the on-shell value a(k,k) — (kcotS — ik)^ 1 where S is 
the elastic scattering phase shift. The solution of Eq. (^) can be obtained 
numerically and is very insensitive to the high momentum modes in the inte- 
gral equation. If, for example, the integral equation is cut off at some finite 
momentum A, the low energy behavior of the solution remains unchanged. 
The EFT in this channel is very predictive. The solution of Eq. (0) is shown 
by the dashed line in Fig. |[ When the corrections of the effective range 
in the two-nucleon interaction are taken into account up to 0((r2/a2) 2 ), one 



obtains the solid curve E3 which 
of van Oers and Seagrave (dots)t 



pees nicely with the phase shift analysis 
I For the scattering length the agreement 

(6.33 ±0.10) 



lVCSt 



,(3/2) 



is even better. The calculation to 0((r 2 /a 2 ) ) gxvt 

fm to be compared with the experimental valuelij a^ 2 ^ — (6.35 ± 0.02) fin. 
Furthermore, the EFT expansion in powers of r^jai converges well, as the 
contributions from 0((r2/a>2) ) to 0((r2/<22) 2 ) to the scattering length are 
a y 3 1 = (5.09 + 0.91 + 0.33) fm, in order. The extension beyond the deuteron 
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Figure 3: fccot5 for fh= 3/2 nd-scattermg to order (r2/a2)° (dashed line) and (r2/a2) 2 
(solid line). The dotslij and the triangleLljj indicate the experimental values. 

breakup threshold and the inclusion of pions in_KSW counting has recently 
been carried out by Bedaque and Griefihammerca 

4 S =1/2 n<i-Scattering and the Triton 

In the 5 = 1/2 channel, the situation is more complicated. The two-nucleon 
interactions can now take place both in the 3 5i and So partial waves. As 
shown in Fig. EL there are two coupled amplitudes, a and b. The amplitude 




Figure 4: Coupled integral equations for S = 1/2 nci-scattering. T (S) dibaryon is indicated 
by grey shaded (black) thick line. 

a which has both an incoming and outgoing dibaryon field T gives the phase 

G 
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shifts for S = 1/2 nd-scattering. However, a is coupled to the amplitude b 
which has an incoming dibaryon S and an outgoing dibaryon T . Although only 
a corresponds to S = 1/2 nd-scattering, both amplitudes have the quantum 
numbers of the triton. From the Lagrangian, Eq. (||), one obtains the coupled 
integral equations for the two amplitudes. After the integration over the time 
component of the loop-momentum and the projection onto the S*-waves has 
been carried out, we have 



!(lA4 + V3p 2 /4- ME) l a{p,k) (9) 

= K{p, k) + - t q \f . K(p, q)[a(q, k) + 3b(q, k)} 
ir Jo Q- k 2 - te 



p2 _ fc2 

2 

i" Jo q 2 - k 2 - ie' 



As before, fc (p) denote the incoming (outgoing) momenta in the center of 
mass frame and ME = 3k 2 / '4 — (l/a 2 ) 2 is the total energy. The kernel 
K(p,q) is given in Eq. (g). a\ (a|) are the scattering lengths in the 3 Si 
^So) NN channels, respectively. The amplitude a(p,k) is normalized such 
that a(k,k) — (kcotS — ik)^ 1 with 8 the elastic scattering phase shift in the 
5 = 1/2 channel. Furthermore, we have introduced a momentum cutoff A in 



the integral equations. Eqs. (|9j, 10) have previously been derived using different 
methodsEH In the limit A — > oo these equations do not have ajinique solution 
because the phase of the asymptotic solution is undetermined^ For a finite A 
this phase is fixed and the solution is unique. However, the equations with a 
cutoff have the same disease as in the boson case: a strmig^ciitoff dependence 
that does not appear in any order perturbation theorvEjtj'Eil The amplitude 
a(p, k — const.) shows a strongly oscillating behavior. Varying the cutoff A 
slightly changes the asymptotic phase by a number of 0(1) and results in large 
changes of the amplitude at the on-shcll point a(p — k). This cutoff depen- 
dence is not created by divergent Feynman diagrams. It is a nonperturbative 
effect and appears although all individual diagrams are UV finite. 

In order to control this strong A dependence, it is useful to note that 
the integral equations ([)], [l0|) are ST/ (4) symmetric in the ultraviolet (UV). 
Therefore it is sufficient to consider the SU(4) limit (a| = a| = (12), since the 
cutoff dependence is a problem rooted in the UV behavior of the amplitudes. 
Furthermore, we note that the equations for a + = [a + b] and a_ = [a — b] 
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decouple in the SU(4) limit. The equations for a + and a_ are given by 



- A (l/aa + ^3p 2 /4 ~ ME) 1 a+{ P ,k) (11) 
n „, 2Jf(A) 2 /" A g 2 dq f . , 2H(A)\ . 



(l/a 2 + ^3p 2 /4 - ME) a_(p,fc) (12) 

2 [ f 
-K(p, k - - / 



g) a-(q,k), 



q 2 ~ k 2 

where we have introduced a contact three-body force H(A) which runs with 
the cutoff A into the equation for <z+. Let us disregard the three-body force 
for a moment. The equation for a + is exactly the same equation as in the 
case of spinless bosons while the equation for a_ is the same equation as in the 
5 = 3/2 channel. The 5 = 3/2 equation is well behaved and its solution is very 
insensitive to the cutoff as discussed in the previous section. Consequently, 
the observed cutoff dependence stems solely from the equation for a + . As 
an example, the cutoff dependence of a + (p,k = 0) is shown by the solid, 
dashed, and dash-dotted curves in Fig. || for three different cutoffs, A = 
1.0, 2.0, 3.0 x 10 4 aa 1 . But we already know the solution, to this problem: a 
three-body force counterterm that runs with the cutoff Alia This is exactly the 
three-body force we have introduced into Eq. (pi]). The dotted, short-dash- 
dotted, and short-dashed curves in Fig. || show the effect of the three-body 
force H(A) on a+{p,k = 0) for A = lO 4 ^ 1 and H = -6.0, -2.5, -1.8. ft 
is clearly seen that the variation of H for a constant A has the same effect 
on the amplitude as varying the cutoff. Consequently, we can compensate the 
changes in the asymptotic phase when A is varied by adjusting the three-body 
force term appropriately. (A more detailed, discussion of the renormalization 
procedure can be found in Bedaque's talkE3) 

We can obtain an approximate expression for the running of if (A) from 
invariance under the renormalization group. Requiring that the equation for 
a + does not change its form when the high momentum modes are integrated 
out, we find 

H(A s _ sin Qo ln(A/A*) - arctg(l/s )) 
1 ' sin( So ln(A/A.) + arctg(l/ So )) ' 1 ' 



where sq w 1.0064E3E3o H(A) contains one new dimensionful parameter, A*, 
which must be determined from experiment. The running of the three-body 
force H(A) according to Eq. ( |l3| ) is shown by the solid line in Fig. ||. The 
dots are obtained by adjusting H(A) such that the low energy solution of Eq. 
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Figure 5: Cutoff dependence of a+(p, k = 0). Solid, dashed, and dash-dotted curves are 
for H = and A = 1.0, 2.0, 3.0 X lO 4 ^ 1 , respectively. Dotted, short-dash-dotted, and 
short-dashed curves show the effect of the three-body force for A = 10 4 a^ 1 and H = 
—6.0, —2.5, —1.8, respectively. 



( |il| ) remains unchanged when A is varied. The observed agreement provides 
a numerical justification for our proceeding. The three-body force is periodic 
with H(A n ) = H(A) for A„ = Aexp(mr/s ) ~ A(22.7) n . Since it enters only 
in the equation for a + = a + b, the three-body force is also SU (4) symmetric. 
Formally, the three-body force term in Eq. (^) is obtained by adding 



2MH(A) ( 

i-3 = — 



A 2 



y g 2 T N\f ■ a)t(f • a)N (14) 
+ \gT9s [N\T -3)\S-t)N + h.c. 
+g 2 s NHS-Ty(S-r)N) , 



to the Lagrangian, Eq. (g). Eq. ( |14[ ) represents a contact three-body force 
written in terms of dibaryon and nucleon fields. Via a Gaussian path integra- 
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Figure 6: Running of H(A) for A* = 0.9 fm 1 : (a) from Eq. (solid line), (b) from 

numerical solution of Eq. Hill) (dots). 



tion it is equivalent to a true three-nucleon force, 



Co = 



2 -^^( K ^(N T r m N)\N^ k a l N)(N T r 2 a l a 2 N) (15) 
1 9t 9 s 



3 2A T 2A £ 



[(N T T 2 a k (T2N^(N^a k T l N)(N T a 2 TiT2N) + h.c 



"4A| 



(N T a 2 TkT2N)\N^TknN)(N T a 2 TiT2N) 



By performing a Fierz rearrangement, it can then be shown that the three 
terms in Eq. (^|) are equivalent. As a consequence, there is only one three- 
body force which is also SU (4) symmetric. Naive power counting would suggest 
that the three-nucleon force scales with l/(Mm*). The three-body force from 
Eqs. ( ]l4| , |l5| ), however, is enhanced by the renormalization group flow by two 
powers of a 2 . Using Eq. (^), it is found to scale as a|/(Mm^) which makes it 
leading order. 

Recently Mehen, Stewart, and WiseEJ found an approximate SU(£) sym- 
metry in the two-nucleon sector for l/a 2 «p« and also noticed that 
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the only S-wave four- nucleoli force that can be written down is St/ (4) sym- 
metric. Furthermore, there are no contact interactions with more than four 
nucleons without derivatives because of the Pauli principle. It is also reason- 
able to assume that the low-energy dynamics of nuclei is dominated by S'-wave 
interactions. Therefore oiir_St/(4) symmetric three-body force together with 
the findings of Mehen et alr3 gives an explanation for the approximate Wigner 
St/ (4) symmetry in nucleiEj 

Now we are in the position to solve the full equations for the broken St/ (4) 
case. Introducing the three-body force from above into Eqs. (Q, |l^), we obtain 



+ V3p 2 /4 - ME) ~* o(p, k) = K{p, k) + (16) 
~ £ q 2 q l k d 2_ lt k (p> k ) + + ^JT^ Hi, k) + b(q, k)\ 



A o 



q J dq 



7r J q 2 — k 2 — 



K(p, q) [3a(q, k) + b(q, k)} + [a(q, k) + b(q, k)} 



ter A_. 

ngthM 



We need one three-body datum to fix the three-nucleon force parameter 
We choose the experimental value for the S — 1/2 nd-scattering le 
4 1/2) = (0.65 ±0.04) fm and find A» = 0.9 fm^ (For the special cutoffs with 
vanishing -ft (A), we recover the results of Refc3.) Although one three-body 
datum is needed as input, the EFT has not lost its predictive power. We can 
still predict (i) the energy dependence of S = 1/2 nti-scattering and (ii) the 
binding energy of the triton. 

The resulting energy dependence of S = 1/2 nd-scattering for three differ- 
ent cutoffs is shown in Fig. ^. It is clearly seen that the introduction of the 
three-body force renders the low-energy amplitude cutoff independent. The 
scattering length is reproduced exactly because it was used to fix A*. The 
agreement for finite momentum is at least encouraging. Our experience from 
the S = 3/2 channel is that the range corrections improve the agreement con- 
siderably (cf. Fig. H) . The dashed curve in Fig. gives a crude estimate of 
these corrections. In the zero range approximation, we have the relation 

VMB 2 = 1/4, (18) 

which holds only approximately in nature. In our calculations we take the 
deuteron binding energy B2 from experiment and determine a| from Eq. 



The dashed curve in Fig. (7] is obtained by taking the experimental value of a\ 
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Figure 7: Energy dependence for 5 = 1/2 nd-scattering with three different cutoffs A = 
1.9, 6.0, 11.6 fm — 1 and A* = 0.9 fm — 1 (hidden under the solid curve). Dashed cupwe gives 
estimate of ranm corrections. Experimental points are from van Oers and Seagrave^j (dots) 
and Dilg et al. Ej (triangle). 

as input and leaving all other parameters unchanged. From the estimated size 
of the range corrections, we anticipate an improved agreement once the range 
corrections are included. (One should also keep in mind that the experimental 
phase shift analysis 113 does not give any error estimate). 

The triton binding energy is obtained from the solution of the homogeneous 
equations corresponding to Eqs. (|l(], [l7]) for E = —B 3 . We find B 3 = 8.0 
MeV for the triton binding energy, to be compared with the experimental 
result B^ xp = 8.48 MeV which is known to very high precision. For a leading 
order calculation, the agreement is very good. The theory without pions seems 
to work for triton physics. However, to draw definite conclusions one has to 
calculate the range corrections for both the energy dependence and the triton 
binding energy. 

In Fig. U we show the bound state spectrum as a function of the cutoff A. 
The shallowest bound state is the triton. Its binding energy is cutoff indepen- 
dent. However, as A is increased new deeper bound states appear whenever 
H(A) goes through a pole. These new bound states appear with infinite bind- 
ing energy directly at the pole. When the cutoff is increased further, their 
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Figure 8: Three-nucleon bound state spectrum for A, = 0.9 fm 1 . The shallowest bound 
state corresponds to the triton. 

binding energy reduces and becomes cutoff independent as well. The poles of 
H(A) can be parametrized as 



One counts n bound states for A n _i < A < A n . However, only the states 
between threshold and A ~ m n ~ l/f2 are within the range of the EFT. By 
solving Eq. ([i"9|) for the number of bound states n, we then obtain 



and recover the well-known Efimov effect E3 In the limit — > oo, an infinite 
number of shallow three-body bound states accumulates at threshold. That 
these bound states are shallow follows from the fact that the binding energy is 
naturally given in units of l/(Afa2i which vanishes as a% — > 00. Furthermore, 
we also recover the Thomas effect til In a hypothetical world where A ~ m, ~ 
I/V2 — > 00, the range of the EFT increases and deeper and deeper physical 



A n = /(A*a 2 )exp(7i7r/s )a 2 

O(l) 



(19) 




(20) 
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Figure 9: Phillips line: S = 1/2 scattering length o,^^ as function of the triton binding 
energy B3: EFT to leading order (solid line); potential models (dots); experiment (cross). 

bound states appear. Consequently, there is an infinitely deep bound state for 
A — * 00. However, for ~ \jv-x as in the real world the deep bound states 
are outside the range of the EFT and their presence does not influence the 
physics of the shallow ones. 

Moreover, the variation of the parameter A» gives a natural explanation 
for the Phillips lineB No matter how good the agreement with the actual ex- 
perimental number is, values of B3 and a^ 2 ^ in different potential models 
are correlated and fall on the Phillips line. Obviously, there is a correlation 
between B3 and 03 . In the EFT framework, different models correspond in 
leading order to different values of A* . Varying A* then generates the observed 
Phillips line. The dynamics of QCD chooses a particular value which up to 
higher order corrections is A* = 0.9 fin - . In Fig. [|we show the Phillips line 
obtained in the EFT compared with results from various potential model cal- 
culations!^ and the experimental values for B3 and a^ 2 ' . Our Phillips line is 
slightly below the one from the potential models. We expect this discrepancy 
to be reduced once range corrections are taken into account. 
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5 Conclusions 



We have studied the three- nucleon system using EFT methods. While for 
nd-scattering in the 5 = 3/2 channel precise predictions are obtained in a 
straightforward way, the S = 1/2 channel is more complicated. It displays a 
strong cutoff dependence even though all individual diagrams are UV finite. 
In this channel a nonperturbative renormalization takes place similar to the 
case of spinless bosons. This renormalization requires an SU (4) symmetric 
three-body force which is enhanced to leading order by the xenormalization 
group flow. Together with the recent results of Mehen et alE-3 this gives an 
explanation for the approximate SU(A) symmetry in nuclei. Furthermore, we 
find that the Phillips line is a consequence of variations in the new dimensionful 
parameter A» which is introduced by the three-body force. A* is not given from 
two-nucleon data alone and has to be determined from a three-body datum. In 
the appropriate limits for the two-body parameters a 2 and we also recover 
the well known Thomas and Efimov effects. 

The leading order of the EFT gives a quantitative description of the triton 
binding energy and the energy dependence for S — 1/2 W-scattering. The the- 
ory shows the potential for a realistic description of the triton once the range 
corrections are included. Immediate applications of the EFT include polar- 
ization observables in nd-scattering and triton properties such as its charge 
form factor. The incorporation of the long range Coulomb force would widen 
the possible applications considerably as pd-scattering and the physics of 3 Hc 
becomes accessible. 



Fina 
counting 



y, since the pionless theory is equivalent to the leading order in KSW 
the success in the three-nucleon system opens the possibility of ap- 



plying the EFT method to a large class of systems with three or more nucleons. 
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